Hyperbolic Models of Homogeneous Two-Fluid 
Mixtures 



S.L. Gavrilyuk and H. Gouin 

Laboratoire de Modelisation en Mecanique et Thermodynamique, Universite 
d'Aix-Marseille, Faculte des Sciences et Techniques, Case 322, Avenue 
Normandie-Niemen, 13397 Marseille Cedex 20, FRANCE 

Yu.V. Perepechko 

United Institute of Geology, Geophysics and Mineralogy, Siberian Branch of the 
Russian Academy of Sciences, 630090 Novosibirsk, RUSSIA 



Abstract. One derives the governing equations and the Rankine - Hugoniot con- 
ditions for a mixture of two miscible fluids using an extended form of Hamilton's 
principle of least action. The Lagrangian is constructed as the difference between 
the kinetic energy and a potential depending on the relative velocity of components. 
To obtain the governing equations and the jump conditions one uses two reference 
frames related with the Lagrangian coordinates of each component. Under some 
hypotheses on flow properties one proves the hyperbolicity of the governing system 
for small relative velocity of phases. 

Sommario. Le equazioni di governo e le condizioni di Rankine-Hugoniot sono 
derivate per una miscela di due fluidi miscibili usando una formulazione estesa del 
principio di minima azione di Hamilton. La Lagrangiana e costruita come differenza 
tra energia cinetica e potenziale dipendente dalla velocita relativa dei componenti. 
Per ottenere le equazioni di governo e le condizioni di salto si usano due sistemi di 
riferimento collegati alle coordinate Lagrangiane di ciascun componente. Nelle stesse 
ipotesi sulle propreta di flusso si prova l'iperbolicita del sistema per piccole velocita 
relative delle fasi. 
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1. Introduction 

The theory of mixtures considers two different kinds of continua: homogeneous mix- 
tures (each component occupies the whole volume of the physical space) and het- 
erogeneous ones (each component occupies only a part of the mixture volume). In 
the second case, the geometrical parameters appear as the unknown variables: the 
volume concentration, sizes of dispersed particles etc. The problem, which is com- 
mon for both types of media, is to describe two- velocity effects, that are responsible 
for the development of instability in mixtures, diffusion of components etc. 

At least three approaches to the construction of two-fluid models are known. The 
most common one for studying of heterogeneous two-phase flows is the averaging 
method (Ishii [1], Nigmatulin [2] and others). Averaged equations of motion are 
obtained by applying an appropriate averaging operator to the balances laws of mass, 
energy etc., which are valid inside each phase. The main problem associated with this 
approach is the closure of the system obtained: the system contains more unknowns 
than equations. Different experimental and theoretical hypothesis are used for the 
closure. Not all of them give well-posed governing equations. For example, it was 
noted by many authors (Drew [3], Kraiko & Sternin [4], Stuhmiller [5] and others) 
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that non-dissipative averaging governing equations of heterogeneous two-velocity 
media are not hyperbolic even for small relative velocity of the mixture components, 
when the equal pressures hypothesis in phases is used. This implies that the Cauchy 
problem for the corresponding equations of motion is ill-posed. The nonhyperbolicity 
was overcome later by many authors: Liapidevskii [6], Ransom & Hicks [7] (two-phase 
separated flow), Biesheuvel & van Wijngaarden [8] (bubbly liquids), Fabre et al [9] 
(slug flow) etc. The well-posedness of governing equations was achieved by using 
additional closure relations for averaged quantities, which are specific to the flow 
type. The hyperbolicity of one-dimensional models only was proved. 

A different approach known as Landau method of conservation laws was ini- 
tially used for constructing models of quantum liquids such as superfluid helium 
(Khalatnikov [10], Landau & Lifshits [11], Putterman [12]). The method consists 
in the following: the requirement of the fullfillment of balance laws of mass, energy 
etc., complemented by the Galilean relativity principle and the Gibbs thermody- 
namic identity fully determines the governing equations of motion. Recently this 
approach was applied to classical fluids (two-velocity hydrodynamics) by Dorovsky 
& Perepechko [13], Roberts & Loper [14], Shugrin [15]. The method does not take 
into account the geometrical characteristics of the mixture components: the vol- 
ume concentrations, sizes of particles etc. In the non-dissipative case it gives also 
hyperbolic models (see, for instance, Khalatnikov [10], where sound velocities for 
superfluid helium are calculated). 

Finally, a third approach called variational method is the most universal. Bedford 
& Drumheller [16], Berdichevsky [17], Geurst [18,19] have applied it for investigation 
of bubbly liquids flows. In particular, Geurst has proven in one- dimensional case the 
hyperbolicity of the governing equations for small relative velocity of phases. 

These methods present three different approaches for description of complex 
media. At present, their common features and distinctions are not quite well under- 
stood. 

We consider the variational approach to describe two-velocity effects in homo- 
geneous mixtures. A physical example of such a flow is a mixture of two miscible 
fluids, or a mixture of two gases with quite different molecular weights. 

In Section 2 we introduce an extended form of Hamilton's principle of least 
action. The Lagrangian of the system is chosen in a general form: it is the difference of 
the kinetic energy of the system, which depends obviously on the choice of a reference 
frame, and a thermodynamic potential, which is a Galilean invariant, conjugated 
to the internal energy with respect to the relative velocity of phases. If it does not 
depend on the relative velocity of components, we have a classical form of Hamilton's 
action for two- velocity systems (see, for instance, the article by Gouin [20], where 
the thermo-capillary mixtures were studied). 

In Section 3 we get from the variational principle the governing equations and 
the Rankine-Hugoniot conditions for shocks. To obtain the desired relations, we used 
two reference frames related with Lagrangian coordinates of each component. 

Conservation laws of the total energy and the total momentum are derived in 
Section 4. We show that, under some restrictions on the flow properties, the gov- 
erning equations admit additional conservation laws in terms of the densities pi, 
P2 and the velocities ui, U2 of components. Without these restrictions the system 
seems not to be conservative. We extend the set of the unknown variables, consid- 
ering the deformation gradients as the required quantities, and rewrite our system 
in a conservative form that gives additional set of possible jump conditions. 

In Section 5, classification of strong discontinuities is done and some difficulties 
of the "right" choice of jump conditions are discussed. 

We investigate in Section 6 the hyperbolicity of the governing system for small 
relative velocity of phases in multi-dimensional case. Under some hypotheses on 
flow properties, we reduce our system to Friedrichs' symmetric form and prove that 
convexity of the internal energy guarantees hyperbolicity of the governing equations. 
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As a convention, in the following we shall use asterisk "*" to denote conjugate 
(or transpose) mappings or covectors (line vectors). For any vector a, b we shall 
use the notation a*b for their scalar product (the line vector is multiplied by the 
column vector) and ab* for their tensor product (the column vector is multiplied 
by the line vector). The product of a mapping A by a vector a will be denoted 
by A < a >. The notation b* A means the covector c* defined by the rule 
c* = (A* < b >)*. The divergence of a linear transformation A is the covector 
divA such that, for any constant vector a 

div Aa. — div (A < a >). 

The letter / will mean the identity transformation, and V will mean the gra- 
dient line operator. The greek indices a, (3 = 1, 2 will denote the parameters of 
components such as the densities p a , the velocities u a etc. 



2. Variational Principle 

Let us suppose that a mixture of two miscible fluids is well described by the velocities 
ui, U2 of two components, the average densities pi, p 2 and the total energy E. 
The total energy is divided into the kinetic energy T and the internal energy U. 
In the following, we will consider only mechanical processes by suppressing thermal 
evolution. Hence, U is purely mechanical part of the total internal energy. The 
kinetic energy, depending on the choice of a reference frame, is represented by the 
classic formula: 

T = g P 1 ! 111 ! 2 + 2 P2 l U2 ! 2 - 

The internal energy U is a Galilean invariant, it does not depend on the reference 
frame. Hence, U depends on pi, p 2 and the relative velocity w = u 2 — Ui. 

Neglecting the dissipative effects, we propose the following extended form of 
Hamilton's principle of least action ( Gavrilyuk et al [21]): 

51 = 0, 1= f'f Ldxdt, L = pi iHil- + p 2 i^J- -W( Pl ,p2,w) (2.1) 
Jt ± Jt> 2 2 

with additional kinematic contraints 

^ + div (piui) = , ^f+div (p 2 u 2 ) = 0. (2.2) 

Here [ti , i 2 ] is a time interval, T) is a fixed bounded domain of the three-dimensional 
space with the boundary dT>. We will suppose that the slipping condition on dT> 
is fulfilled for each component. 

The internal energy U is the partial Legendre transformation of the potential 
VF(pi,p2,w) with respect to the variable w: 

U = W(p 1 ,p2,w) -— W=VV+1 W, 1 = , 

0W ow 

BU 

W = U-i*w, w* = (2.3) 
oi 

The vector variable i, which is also a Galilean invariant, can be called the relative 
momentum. Below, we will suppose that U is invariant under rotations (the case 
of isotropic media). Hence, U is a function of i = \ i | and W is a function of 
w = | w |. In this case the relations (2.3) can be represented in the form, 

U = W - — — w = W + i to, i = — — , W = U - i w, w = —. (2.4) 
out OW Ol 
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If U does not depend on i, the variational principle (2.1), (2.2) coincides with 
classical Hamilton's principle of least action: the Lagrangian is the difference of the 
kinetic and the internal energy. The relation (2.3), (2.4) between U and W will 
be justified later when the conservation law of the total energy will be obtained. 



3. Governing Equations and Jump Conditions 

Let x be Eulerian coordinates, X a be Lagrangian coordinates of the a-th com- 
ponent, a = 1,2. The relation between Lagrangian and Eulerian coordinates is 
given by the diffeomorphisms of the domain D into T>: 

x = <P a (*,X„), X Q = ip a (t,x), (j> a o ip a = I, a = 1,2. (3.1) 

Let 

F a = j&-, a = l,2 ( 3 - 2 ) 

be the deformation gradient at X a (or Jacobian matrix of the mapping <j> a ). Let 
us define the virtual motions of the mixture (Serrin [22], Gouin [23] and others): 

x = * a (t, X a , e a ), X a = * a (t, x, e a ), * Q o * a = J, (3.3) 

where e a varies in the neighbourhood of zero. The real motion corresponds to 
Ea = 0: 

* Q (t, X Q ,0) = 4> a (t, X a ), * a (t, X, 0) = j> a (t, X). 

The associated variations <5 a x and S~K a are defined by the relations: 

8 a x = 1*2 (t,X a ,0) , <5X Q = (t, x , ). (3.4) 

The definitions (3.1) - (3.4) imply that 

5 a x = -F a < S X a > . (3.5) 

Let 

/a(*,x), / a (t,X a ) = / a (t,0 a (t,X a ) ), 

A ~ A 

/a (t,X,£ a ), / a (<,X a ,£ a ) = / Q (<,* a (t,X ai £„),£„) 

be the unknown quantities of the a-th component (such as the density p a , the 
velocity u Q etc.) in Eulerian and Lagrangian coordinates and their perturbations, 
respectively. One defines Eulerian and Lagrangian variations of the variable f a : 

A ~ 

d f ° d f 

5 f" = 7T^ (*- x >°) - 3 f« = (*,X a ,0). 

O £a O £a 

It yields (Berdichevsky [17], Gouin [20,23]) 

° Bf 
S f a = 5 f a + -JZ SaX . (3.6) 

Using the Euler formula 

' d 8 a x 



<5 (det F a ) = det F a div(S a x) , div(8 a x) = tr 
and the mass balance (2.2) in the form 

°p a det Fa =°p a (0,X a ), 



8 x 
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we obtain Lagrangian variations of p a : 

8 °p a = - °p a div S a x. (3.7) 



We also note that 



S u a = A S a x. (3.8) 



Using (3.6), (3.7), (3.8), we get 

Sp a = - div (p a c5 Q x) , 8 u a = ^- S a x - < <5 Q x >, (3.9) 

d f) 

where V* means the gradient column operator. We assume here that 8 a n are 
functions of Eulerian variables. We define the vectors 

* 1 dL . 1 dW w* 
K a = — - — = u a - (-1) — , (3.10) 

pa OUa pa OVJ W 

where 

1 2 1 2 

L(pi, p2, Ul, U2) = - pi |ui| + - p'2 |U2 1 — (pi,p2,w) ,W = I U2 — Ul|. 

Varying the variables (pi,ui) and (p 2 ,u 2 ) independently and denoting by 5 a T 
the corresponding variation of the functional T , we obtain 

= I j j dp,, ( i|u Q | 2 - ) + K, (1'u,, j r/x r// 



ti jo 



Taking now into account the formulae (3.9), we get 
S a I = I I — di« (p a <5 a x 



1 , ,2 

- U a - 

2 Opa 



+ p a K* ^ ^ <S a x - < <S a x > J J dx di = 

+ V [dp-a — J J + «(^ K «'- X ) + 

+ div ^p a u a K; + p a (j!?- - < <S a x > ^ j dx di = 0. (3.11) 

The mapping I means here the unit tensor and, as previously defined in the 
Introduction, for any vectors a and b, a* b is the scalar product and a b* 
is the tensor product. If all the functions in (3.11) are smooth in the domain T> 
and the variations 5 a x vanish on dT>, the divergence terms do not play any role. 
Therefore, we obtain the equations of motion: 

dK a 9K a <9u* (dW |u„ 

^r + ^r <Ua>+ ^r <Ka>+v [dp- ~ — 



= 0, (3.12) 
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or, equivalcntly 
dK 



dt 



+ rot K a x u a + V* - + K* u 



aw 



(3.12') 



If the function W (pi, pi, w) (or the internal energy U(pi, p2, i)) is given, 
the equations (3.10), (3.12) with the mass conservation laws (2.2) form the closed 
system of the governing equations. We will show in Section 6 that, under natural 
restrictions on the internal energy U (or W), the system is hyperbolic if the relative 
velocity is sufficiently small. 

Now, suppose that the domain T> x [ t\,t2 ] is divided by a singular surface 
S(t) having at any of its points the normal unit vector n and the normal speed of 
displacement D n . Suppose also that at any point of S(t) the right and left limits 
of K„ and p a exist, but not necessary equal. Then, the divergence terms in (3.11) 
give the jump conditions (the Rankine-Hugoniot conditions): 



K* + p a [ q — 



dW 
dp 



D n p a K* <5 a x 



< <5 Q x > - 







(3.13) 



where the square brackets denote the jump. 

If the singular surface S(t) is a Shockwave for a-th constituent (i.e. n* u a — 
D n ^ which means that the particles cross the surface), the formula (3.13) can 
be symplified. Indeed, using, (3.5) and the fact that [ <5X a ] = 0, we obtain: 



Pa ( n* U a - D n )K a F a + Pc 



( dW 
\dp a 



* 7-1 

n F a 



0. 



(3.14) 



The equations (3.14) contain not only ~K a and p a but also the deformation 
gradient F a . However, in the next section, we will obtain shock conditions in terms 
of K. a and p a , using a conservative form of the governing equations in Eulerian 
coordinates. 

Finally note that the governing equations could be also obtained by using the 
method of the Lagrange multipliers. 



4. Conservation Laws 

Conservation laws, i.e. the expressions of the form 

d Pq 



dt 



+ div P = 0, 



where Po, P are functions of unknown variables, play an important role in the 
theory of hyperbolic equations (see, for instance, text-books by Serre [24] or Smoller 
[25]). The property of conservativeness of mathematical models, when the number 
of linear independent conservation laws admitted by the model is not less than the 
number of unknown variables, is necessary to determine weak solutions of the system. 
Some conservation laws can play a role of entropy, i.e. all admissible solutions of the 
system of conservation laws must satisfy the " entropy inequality" 

^ + div h < , 
a t 

where ho, h are functions of the unknown quantities. 
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The equations of motion (2.2), (3.12) admit two obvious additional conservation 
laws of the total momentum and the total energy, corresponding to the invariance 
of the Lagrangian 

1 2 1 2 

L = g P 1 l Ul + 2 P2 |u 2 | - W (pi , P2 ,w) 

with respect to time and space shifts. They can be obtained either using the theorem 
of E.Noether (Olver [26], Ovsyannikov [27]) or by direct calculations. The momen- 
tum conservation law is obtained multiplying the equations (3.12) by p a and then 
summing: 

d . J . / 

— ( pi Ui + p2 U 2 ) + dlV I pi Ui Uj + p 2 u 2 u 2 - 

aw ww* / aw aw T \ r \ „ 

- 75 + Pi o— + 75 W 7 = 0. 4.1 

aw w y dpi dpi ) ) 

It admits also an alternative form in terms of K Q (see the definition (3.10)): 

A ( Pl ki + P2 k;) + 

+ dl« I pi Ul Kt +p 2 U 2 K 2 + I pi — + p 2 -7^ ~ Wl J 1 = 0. (4.1') 

Multiplying (3.12) by p a u a and then summing, we obtain the energy conser- 
vation law 

d ( |ui| 2 |u 2 | 2 „, dW\ J . ( ( |ui| 2 dW 

/ |u 2 | 2 dW \ dW ( „ „ \ w \ „ IA „, 

+ ^ Ul ( v ^ + aplJ-^7( v U2U2 " UlUl J < ™ >J =°" (42) 

It admits an alternative form: 
d ( |ui| 2 |u 2 | 2 rr \ / /aw |ui| 2 ^ \ 

a7 ("V" + " 2 Sr + + (/ lUl ( apT " + KlUl ) + 

+ P2U 2 (^-M + K|U2))^0. (4.2') 

The energy conservation laws (4.2), (4.2') explains now the relation (2.4) between 
U and W: 

aw 

U = W-w 

aw 

It follows from (3.12'), that an additional conservation law can be obtained, if 

rot K a = 0. (4.3) 

In that case 

§? + v(«- + £-¥)-' (4 „ 

The equations (4.3), (4.4) are compatible. Indeed, the equation (4.4) admits a con- 
sequence 

— rot K a = 0, 

a t 
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and hence, (4.3) can be considered as a restriction for the initial data. 

To our knowledge, there is no additional conservation laws in terms of the vari- 
ables p a and u Q . Finally the system (2.2), (3.12) contains eight desired variables 
p a , u a , ct = 1,2. If rot K. a 7^ 0, it admits only six conservation laws (2.2), (4.1), 
(4.2). Hence, in the general case, the system seems not to be conservative. 

Below, we will extend the set of desired variables, considering the deformation 
gradients F a as the unknown quantities. We will show that the extended system 
is a system of conservation laws. 

Straightforward calculations show that F a satisfies the equation: 



Using the Euler formula 

% ( det F a ) = det F a div u a . % = ^- + u a V , 
at dt at 

we obtain 

d a ( F a \ ( 8u a 



\ i , , i - \ ■ div u a I ) ■ (4.6) 
dt \ det F a J \ ox J det F a 

Taking into account (4.5), we rewrite (4.6) in the divergence form 

d ( (Fa)T \ + _d_ ( (ug) k (Fg)T - (UoT (F a ) k s \ =Q (4 ?) 



dt \ det F a ) dx k \ det F a 

Here (F a )T are the components of F a (m denotes the lines and s denotes 
the columns), (u„)' are the components of u a - The repeated latin indices imply 
summation. Conversely, applying the operator div to the equation (4.7), we obtain: 



^_ div , = 0. 

dt \ det Fa J 

Hence, we can replace the equation (4.5) by the evolution equation (4.7), considering 
(4.5) as the restriction for initial data: if the condition (4.5) is fullfilled at t = 0, it 
is valid for any time. We note that the divergence form (4.7) was earlier obtained 
by Godunov and Romensky [28] in the theory of elasticity. 

Finally, by using (4.7), we get the equation (3.12) in the conservative form: 

9 fK « F ^ + ^((uaKa + (^-^)l)^)=0, (4.8) 



dt \ det Fa J \\ \dp a 2 J J det F a 

which represents the conservation of local momentum of a-th phase. In the next sec- 
tion, we will show that the conservation law (4.8) corresponds to the jump conditions 
(3.13) obtained from the variational principle (2.1). 

The equations (2.2), (4.5), (4.7), (4.8) are in the conservative form relative to 
the variables p a , u a and F a - It is also shown that the equations admit also 
conservation of total momentum, conservation of total energy and , in the case of 
rot ~K a — , conservation of K Q . 



5. Analysis of Rankine-Hugoniot Conditions 

Let S(t) be a singular surface with the unit normal vector n and the normal 
velocity D n , where the functions p a , K Q and F a have jumps. The equations 
(2.2), (4.7), (4.8) imply the following Rankine-Hugionot conditions: 

[ P a ( n* u a - D n ) ] = , (5.1) 
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( n* u Q - D n ) K 



[ ( n* u a - D n 



det F a 



det F a 
dW 

dpa 



det F ( 

!,, |2 



o , 

Fa 



det F a 



9 

(5.2) 
= 0. (5.3) 



If n* u Q — D n — 0, we call it a contact discontinuity. If n* u„ — D n 0, we 
call it a shock wave. 

Let us consider the case of contact discontinuity. It follows from (5.2) - (5.3) that 



dW 

dpa 



u Q n* F a 

det Fa 

luJ 2 



= 0, 

n* F a 
det F a 



(5.2') 



0. 



(5.3') 



Multiplying (5.2') by the normal vector n* and taking into account that [n*u a ] 
0, we get 

" n* Fa 



det F a 



0. 



(5.2") 



The equations (5.2') and (5.2") imply that the velocity u Q is continuous. Hence, it 
follows from (5.2'), (5.3') that the conditions on the contact discontinuity are given 
by 

dW 



[ U a ] = 0, 



dpa 



= 0. 



(5.4) 



Let us consider shock waves. First of all, we note that the condition (5.1) can be 
rewritten in the form 

[ p a det Fa ] = 0. (5.1') 

Hence, the relations (3.14) derived from the variational principle coincide with (5.3). 
As in the case of the contact discontinuity, we rewrite the equations (5.1) - (5.3) 
in terms of u a , pa- Multiplying (5.2) by n*, we get the equation (5.2"). Further, 
let S(t) be a shock surface in Eulerian coordinates, and S a (t) be its image in 
Lagrangian coordinates of the ct — th component. Let q a be a tangent vector to 
S a (t). Then, q = F a < q Q > is a tangent vector to S(t). Multiplying (5.3) from 
the right by q a and taking into account that n* q = 0, we obtain: 



(n* u a - D n ) 



K* q 

det F a 



0. 



It follows then from the last relation and from equations (5.1), (5.1') that the tan- 
gential component K. aq of K Q = n (n* K a ) + K aq: U K a q — 0, is 
continuous: 



K 



0. 



Multiplying (5.2) by K* g , we have: 



K* p 
aq r ol 



(n* u Q 



D n 



det F a 



(KaqUa) U* F a 



det F a 



(5.5) 
(5.6) 



Replacing in (5.3) K Q by K Q9 + n (n* K a ) , we get 



(n* u a - D n ) 



+ 



dW 

dpa 



(K* 9 + n* (n* Kg) ) F a 

det F a 

n* F a 



+ 



det F a 



0. 



(5.7) 
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Equations (5.2"), (5.6), (5.7) imply then that 
dW | u, 12 

dp a : 



+ K* u Q 



0. 



(5.8) 



Equations (5.1), (5.5), (5.8) are the Rankine-Hugionot conditions for shocks in terms 
of the variables u a ,p a . It is worth to note that the jump conditions (5.5), (5.8) for 
shocks coincide with the jump conditions for the equation (4.4). Nevertheless, we 
did not use in our derivation the hypothesis rot K Q = 0. 

The conservation laws (4.1), (4.2) imply also additional jump conditions for the 
total momentum and energy. Hence, we obtain an overdetermined system of the 
Rankine-Hugoniot conditions in terms of the physical variables p a , u a , a = 1,2. 
This is a consequence of the well-known fact that the same system of equations can 
be written in different divergence forms each of which defines different weak solu- 
tion (see, e.g., [24], [25]). We can now question, which divergence form is the more 
appropriate? For one-velocity systems this choice is unambiguous. For example, for 
isentropic gas flows we use the conservation of mass and momentum. The mechanical 
energy plays the role of entropy: it decreases through the shocks [24]. The choice of 
appropriate shock conditions for the two- velocity case is less clear. Hamilton's prin- 
ciple provides a set of Rankine-Hugoniot conditions (3.14). As was shown in Section 
4, these last conditions correspond to the divergence form (4.8) which represents the 
conservation of local momentum. Formally, equations (5.5), (5.8), which are issued 
from (3.14) and supplemented by the equations of mass conservation (5.1), form a 
complete set of Rankine-Hugoniot conditions. Similar to the one-velocity isentropic 
case, the energy conservation law (4.2) (or (4.2')) should apparently play for shocks 
the role of "entropy" inequality: 



-D n [ E ] + 



E 



f dW 

\ dpa 



+ K* u Q 



< 0, 



Pi J 



ui 



+ P2- 



U2 



+ u. 



2 r 2 

The jump conditions obtained are inconsistent with the conservation of the total 
momentum (4.1) (or (4.1')): 



2 

E 

a=l 



p a (n* u a - D n ) K* + n* 



Pi 



dW 
dpi 



+ P2 



dW 
dp2 



W 



= 0. 



Finally, we note that the system of the jump conditions for two-fluid models is 
generally underdetermined. This does not permit to define weak solutions. In our 
case, this system is overdetermined. Hamilton's principle provides a complete set of 
appropriate jump conditions. Is that choice correct? Only physical arguments are 
able to give a definite answer to this question. 



6. Hyperbolicity 

The property of hyperbolicity of governing equations is very important, because it 
implies the well-posedness of the Cauchy problem. Below, we will give a sufficient 
condition of the hyperbolicity of the system (2.2), (3.12) in the multi-dimensional 
case provided that rot K a = 0. 

First, we transform our system to a symmetric form. Considering the Lagrangian 
of our system 

2 

l = E 2 Pa ' Ua ' 2 ~ w ,P2 

a=l 



Meccanical998.TeX.tex; 11/02/2008; 1:13; no v. ; p. 10 



Hyperbolic Models of Homogeneous Two-Fluid Mixtures 



11 



we get 

2 



i dL 

^- Pa du a ] = 



, T / dL , dL , \ ( dL , ] 
dL= } - — dp a + ^ — du a = > ■= — dp a H 

a — l x ' a — 1 

Q = l ^ ' 

= ( ( - K *a J d p a + K* a d(p a U a ) J = ^ (tT Q dp a + K* dja ) , 



where 



dL 



xr* ( ' 

- K„ u a = - - 



dp a \dp a 2 

Or, equivalently, 

2 2 

d (L — cr a p a ) = ^ - p a da a + K* dj a . (6.1) 

(x — l a— 1 

Let us introduce 

2 2 ^ 

G((7i, (72, jl, j 2 ) = P2, jl, J2) - ^2 a a p a = L - ^ p Q . (6.2) 

a—l a—l 

The function G is a partial Legendre transformation of L (pi, p 2 , ji, $2) with 
respect to the variables p a : 



dG dG . 

a = - P«> "S^ - = K a (6.3) 

da a dj a 



By using (6.1) - (6.3) we get: 



d f dG 
dt \ da a 

d_( dG_ 

dt \ d'ia 



div ja = 0, 
V a a =0. 



Or 



dt \ da a I V 9a 



2 



0=1 

2 



i(£)-*(£(s>"JJ-' 1 

The system (6.4), (6.5) can be rewritten in a symmetric form (Friedrichs [29], 
Friedrichs & Lax [30], Godunov [31], Godunov & Romensky [28]): 

A W + B% £ = ° ' A = A * ' B ' ={Br ' i = 1 > 2 ' 3 ' ( 6 - 6 ) 
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where 

f\ 2 r 

U = (<7i , <7 2 ,Ji, J 2 ) , A = — T 

and the matrices £T can be obtained from (6.4), (6.5). If, moreover, the matrix A 
is positive, the system (6.6) is hyperbolic. It is worth to note that equations (6.4), 
(6.5) admit a conservation law of the form 

|( E ^§- a + f>- G ) - dra ( E - °- 

which coincides with the equation of energy (4.2'). (see also Godunov [31]). 

So, we need to prove convexity of G (ai, <t 2 , ji, J2). As we have previously 
mentioned, G is the Legendre transformation of L (pi, p2, ji, J2 ) with respect 
to pi, p2 (see the formulae (6.2) , (6.3)). If L is a convex function with respect 
to ji, J2 and concave with respect to pi, p 2 , then G (<7i, 02, ji, j2 ) will be 
convex, that means the hyperbolicity of our system. Hence, it is sufficient to prove 
that the symmetric matrices Ly 3 and L pp , defined below, are positive and negative 
definite, respectively: 

/jj I 0. 



/ Aii.ii 


Ahj2 


V Aiij2 




f L P1P1 


LpiP2 




LpiP2 



Lpp = I I < 0. 



It follows from (1.2) that for the isotropic case 

L{puP2, jija) = ]{ ^+ ] ^-W[ l n. l>1 . 
2pi 2p 2 



J2 _ Jl_ 

P2 Pi 



Then, 



T _ jc* , n « 1 9W w* 

^Ja — — ^ — -L J "5 ! 

pa pa OW W 

1 d 2 W WW* 

Li 



J1J2 „„ fl,„2 



plp2 9to 2 
L 1 1 9 2 iyww' 

JqJq p a p 2 <9w 2 

Straightforward calculations shows that Lii > 0, if — — tt < 0. 

dw 2 

Now, we calculate L pp : 

0L_ = _ _ 8W + W W*jg 

dpa 2p a 2 dp a dw W pi ' 

dW w*ji 

9to wpf 

w*j Q 



d 2 L 


d 2 W 


- —1 


dpidp2 


dp\dp2 


dp2 \ 


d 2 L 


Li. | 2 d 2 w 


d 


dpi 


pi dpi 


dp a 



dw wpi 

Consequently, L pp < if the velocities u Q are sufficiently small, and the 
function W is convex with respect to pi, p 2 . However, the governing equations 
are invariant under the Galilean group of transformations 

x' = X + Ut, U a = U a + U, t' — t. 
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That means that the condition "the velocities u„ are sufficiently small" can be 
replaced by "the relative velocity w is sufficiently small". Hence, the conditions 

d 2 W Q d 2 W Q d 2 W d 2 W _ ( d 2 W \ 2 Q 
dw 2 ' dp 2 ' dpi dp\ \dp\dp2 J 

guarantee the hyperbolicity of our system for small relative velocity of phases. Due 
to (2.4), the inequalities (6.7) mean the convexity of the internal energy U(pi, P2,i), 
that corresponds to a natural condition of thermodynamic stability. 

Finally, we have established that the thermodynamic stability implies the hyper- 
bolicity of the governing equations for small relative velocity w, provided that 
rot Ki = 0. The last condition is always fullfilled for one-dimensional flows. 
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